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Abstract 

In  this  study,  we  explored  the  possibility  of  enhanced  spontaneous  emission  of  radiation  beyond 
the  free  space  value  by  analyzing  a  semiconductor  superlattice  structure  placed  in  a  microcavity 
whose  resonant  modes  were  tuned  to  the  Bloch  frequency.  In  particular,  we  considered  the  spon¬ 
taneous  emission  of  Bloch  radiation  into  the  dominant  mode  of  a  rectangular  waveguide.  In  the 
analysis,  the  quantum  radiation  field  was  described  by  the  waveguide  quantized  electromagnetic 
field  in  the  Coulomb  gauge,  and  the  instantaneous  eigenstates  of  the  Bloch  Hamiltonian  were  used 
as  basis  states  to  analyze  the  Bloch  dynamics  to  all  orders  in  the  constant  external  electric  field. 
The  results  predict  that  the  spontaneous  emission  occurs  with  frequencies  equal  to  integral  mul¬ 
tiples  of  the  Bloch  frequency  without  any  ad  hoc  assumptions  made  concerning  the  existence  of 
Wannier-Stark  ladder  levels;  such  quantization  effects  arise  from  a  natural  consequence  of  the  im¬ 
plicit  quantum  selection  rules.  It  was  shown  that  the  power  radiated  into  the  dominant  mode  of  a 
rectangular  waveguide  can  be  enhanced  by  an  order  of  magnitude  in  comparison  with  that  for  the 
free  space  spontaneous  emission  by  tuning  the  Bloch  frequency  to  align  with  the  spectral  region 
of  the  waveguide  spectral  density  peak.  For  GaAs-based  superlattices,  the  power  radiated  from 
spontaneous  emission  due  to  Bloch  oscillations  in  the  terahertz  frequency  range  was  estimated  to 
be  about  several  microwatts. 
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I.  Introduction 

Electrons  in  a  miniband  of  a  semiconductor  superlattice  (SL),  accelerated  by  an  electric 
field  E  applied  in  the  transverse  direction  to  the  SL  layers  without  scattering,  can  undergo 
Bloch  oscillations1-3  with  frequency  ujb  —  eEa/h ,  where  a  is  the  SL  period.  The  Bloch 
frequency  ujb  can  be  tuned  with  external  field  E  so  as  to  vary  from  the  high  gigahertz 
to  the  terahertz  (THz)  frequency  range,  thus  suggesting  the  applicability  of  such  internal 
oscillations  as  a  tunable  solid-state  source  for  submillimeter  wave  radiation.2  Most  recently, 
Bloch  oscillations  have  been  confirmed  by  using  time-domain  terahertz  emission  spectroscopy 
in  GaAs-based  SL  structures.4,5  There  have  been  also  a  lot  of  discussions  relating  the  coupling 
of  coherent  Bloch  oscillations  to  other  fundamental  solid-state  excitations  such  as  photons,6’' 
optical  phonons,8  and  plasmons.9  In  this  study,  we  analyze  the  spontaneous  emission  (SE) 
of  radiation  resulting  from  a  Bloch  electron  accelerated  by  a  constant  electric  field  through 
a  miniband  of  a  SL  structure  placed  in  a  resonance  microcavity. 

The  use  of  an  electromagnetic  cavity  is  well  known  to  maintain  resonant  conditions  at 
a  given  frequency  to  provide  a  positive  suitable  feedback  in  various  microwave  oscillator 
schemes10  as  well  as  in  optical11  and  terahertz12  laser  systems.  An  increase  of  the  overall 
emitted  THz  power  of  more  than  one  order  of  magnitude  has  been  reported  by  placing  a 
surface- field  emitter  inside  a  THz  cavity.13  As  axamples  relevant  to  this  work,  many  efforts 
have  been  focused  on  increasing  of  the  SE  rate  in  optical  microcavities.14  The  effect  of 
electromagnetic  wave  confinement  (i.e. ,  when  at  least  one  dimension  of  the  cavity  is  of  the 
order  of  the  radiative  wavelength)  is  to  redistribute  the  free-space  mode  spectral  density  so 
as  to  increase  it  at  some  frequencies  and  to  decrease  it  at  others.  Therefore,  for  an  active 
medium  in  a  cavity  structure,  the  SE  rate  can  be  enhanced  or  diminished  depending  on  the 
position  of  the  emission  frequency  relative  to  the  cavity-mode  spectral  density.14,15 


II.  Bloch  oscillator  model  configuration 

We  assume  that  the  SL  structure  is  placed  into  a  waveguide  with  rectangular  cross  section 
Lx  x  Ly  and  length  Lz,  where  the  coordinate  axes  are  chosen  to  be  along  the  waveguide  edges. 
The  dc  electric  field  E  is  applied  along  the  y  axis,  which  is  also  the  SL  growth  direction. 
The  electromagnetic  field  inside  the  waveguide  with  assumed  perfectly  conducting  walls  is 
determined  by  the  guided  modes  corresponding  to  standing  waves  with  respect  to  the  X 
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and  Y  axes  [designated  by  an  integer  pair  (m,  n)]  and  propagating  waves  along  the  Z  axis 
characterized  by  propropagation  constant  qz.  Such  modes  form  a  complete  and  orthogonal 
basis  set  for  describing  the  electromagnetic  fields  within  the  waveguide.  In  the  following,  we 
will  consider  only  transverse  electric  (TE)  modes,  where  the  electric  field  is  perpendicular 
to  the  direction  of  propagation.  For  practical  cases,16  the  most  important  of  all  confined 
modes  is  the  TEio  mode  (m  =  1 ,  n  =  0),  which  is  the  dominant  mode  of  a  waveguide 
with  Lx  >  Ly.  This  mode  has  the  lowest  attenuation  and  its  electric  held,  Er,  for  the 
chosen  system  geometry,  is  polarized  in  the  direction  of  the  dc  held  E.  We  consider  only 
one  excited,  TE10,  mode  while  all  the  other  less  effective  TE  and  TM  modes  are  ignored. 


For  this  mode,  we  have  the  following  zero  electric  held  Er  x  =  Er,z  =  0  and  magnetic  held 
Hr,y  =  0  components.  Thus,  the  vector  potential  Ar  of  the  held  has  only  one  nonzero  (y) 
component,  so  that 


Ar,y  — 


I  Anhc2 
ujqeV 


sin(qxx)  ( aqeiqzZ  +  a)qe  lQzZ)  , 


where  qz  =  2nl/ Lz  (l  is  an  integer),  qx  =  n/Lx,  c  is  the  velocity  of  light  in  vacuum,  V  = 
LxLyLz  is  the  waveguide  volume,  ceq  and  aq  are  the  creation  and  annihilation  boson  operators, 
and  £  is  the  dielectric  constant  of  the  medium  filling  the  waveguide.  The  normalization 
constant  in  Eq.  (1)  is  chosen  in  such  a  way  that  the  Hamiltonian  for  the  quantized  radiation 
held  has  the  usual  form  Hr  =  J2qhtuqal]aq,  where  uq  =  cnc[l  +  ( qz/qx )2]1//2  is  the  mode 
dispersion  relation,  and  uc  =  qxc/ \fe  is  the  angular  cutoff  frequency.  The  guided  mode 
wavelength  is  written  as  A  =  Ac/^u^/ciy)2  —  l]1/2,  where  Ac  =  2 Lx  is  the  cutoff  wavelength.10 


III.  Bloch  Hamiltonian  and  quantum  dynamics  based  on  instantaneous  eigen¬ 
states 

In  this  analysis,  the  dynamical  properties  are  considered  for  the  situation  in  which  the 
electron  is  confined  to  a  single  miniband  ‘no”  of  a  SL  while  the  effects  of  interband  coupling17 
and  electron  intraband  scattering  are  ignored.  Therefore,  the  quantum  dynamics  is  described 
by  the  time-dependent  Schrodinger  equation 

«^l*n0(«))  =  ff|S'„„(i)>.  (2) 

where  the  exact  Hamiltoian  H  can  be  reduced  to  a  sum  of  the  following  Hamiltonians 
H  =  Hq  +  Hr  +  Hj.18  Here  the  first  two  terms  represent  the  Hamiltonian,  H0(t)  = 
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[p  +  Pc(£)]2/2mo  +  Vc{t),  for  a  single  electron  in  a  periodic  crystal  potential  Vc{r)  inter¬ 
acting  with  a  homogeneous  electric  field,  and  the  Hamiltonian,  Hr,  for  the  cavity  mode 
electromagnetic  field.  The  Hamiltonian  Hj(t)  =  — (e/m0c) Ar  ■  [p  +  pc(t)],  for  the  first-order 
interaction  between  the  quantum  field  and  the  Bloch  electron,  couples  both  subsystems  H0 
and  Hr,  and  causes  transitions  between  the  accelerated  Bloch  electron  states  through  photon 
absorption  and  emission;  pc(f)  =  e  E (t')dt',  mo  is  the  free  electron  mass.  The  total  vec¬ 
tor  potential  consists  of  A  =  Ac  +  Ar,  where  Ac  and  Ar  describe  the  external  electric  field 
and  the  cavity  mode  quantized  radiation  field,  respectively.  Below,  starting  with  the  total 
Hamiltonian  H,  use  is  made  of  first-order  time-dependent  perturbation  theory  to  calculate 
SE  transitions  probabilities  between  states  of  H0  +  Hr  while  regarding  Hj(t )  Ar-[p+pc(t)] 
as  a  perturbation.19  The  solution  to  |\I/no(t))  of  Eq.  (2)  can  be  represented  in  terms  of  eigen¬ 
states  of  basis  states  |^nokp),  {riq,s})  =  |VVik(t))  |{'^q,s})  °f  the  unperturbed  Hamiltonian 
Hq  “l-  Hr  clS 


l^no  «>  =  EE  v,.  s}  (^)  t)  |t/;nok(i) )  {'^q,s} ) 
Xexp  jf  [£no(k(t'))  +  ^  ^9nq!s]^,|  > 


(3) 


where  the  summation  over  k  is  carried  out  over  the  entire  Brillouin  zone,  and  {nq>s}  is 
specified  over  all  possible  combinations  of  photon  occupation  number  nChS  with  photon 
wave  vectors  q  and  polarization  fqjS.  The  instantaneous  eigenstates  of  H0  are  given  by1' 
t/^n0k(i)  (r,  t)  =  f2”1/2e*K'ru„ok(q(r,  t) ,  where  Mn0k(t)(r,£)  is  the  periodic  part  of  the  Bloch  func¬ 
tion,  k (t)  =  K  +  p c(t)/h,  and  the  values  of  the  electron  wave  vector  K  are  determined  by 
the  periodic  boundary  conditions  of  the  periodic  crystal  of  volume 

For  the  one-photon  SE,  which  assumes  that  initially  no  photons  are  present,  the  proba¬ 
bility  amplitude  in  the  wave  function  of  Eq.  (3)  satisfies  the  initial  condition  A{n  j(k,  to)  = 
{5nq  sio}  Ko  at  time  t  —  t0  when  the  electric  field  is  turned  on.  Here,  K0  and  nq  =  0 
are  the  initial  values  of  K  and  nqiS.  The  probability  amplitude  for  SE,  A9(k0,t),  is  now 
evaluated  in  first-order  perturbation  theory  as 


Aq(k0,t)  =  D(qx/q)1/2  /  dt' vy(k0  -  q,) 


1 10 


xexp 


h 


£no(k0)  -  £no (k0  -  q^)  -  tiuq]dti  \  , 


't0 


(4) 
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where  D  =  —i^/nca/u3c£V,  a  =  e2 /he  is  the  fine  structure  constant,  k0(t)  =  K0  +  pc(t)/K, 
q j  =  {±qx,  0,  qz}  with  ”  +  ”  for  j  —  1  and  ”  —  ”  for  j  =  2,  and  q  =  (g2  +  g2 )1/2.  Then  the 
emission  process  results  in  the  well-known  SE  probability 

(5) 

q  3= 1,2 


IV.  Quantum  selection  rule 

In  evaluating  Aq(k0lt),  we  take  into  account  that  the  external  dc  field,  E.  is  along  the  Y 
axis;  then,  it  follows  that  k0y(t)  =  K0y  +  eE(t  —  t0)/h.  In  taking  advantage  of  the  periodic 
properties  of  the  terms  in  Eq.  (4),  x4g(k0,t)  is  evaluated  in  clocked  integral  multiples  of  the 
Bloch  period,  so  that  t  =  Ntb,  where  tb  =  27 t/lub,  the  time  to  traverse  one  period  of 
the  Brillouin  zone.  The  integral  in  Eq.  (4)  over  time  can  be  replaced  by  an  integral  over 
koy  through  the  substitution  dt  =  ( U/eE)dkoy .  Then  the  probability  amplitude,  at  integral 
multiples  of  the  Bloch  period,  can  be  expressed  through  that  over  the  single  Bloch  period, 
Tg.1('18  Thus  we  obtain 


K(ko,iVTB)|2  =  l^(ko.rB)|2, 

sm  yPq/ Z) 

where  the  parameter  f3q  is  given  by 

fjj  \  rKOy+Gy 

pq  =  2vr  —  +  —p  dk0y  [eno(k0)  -  £no(ko  -  qy)]  , 

Vb  eE  JKoy 


(6) 

(7) 


and  Gy  =  2tt /a,  the  y  component  of  the  SL  reciprocal- lattice  vector. 

From  Eq.  (6),  it  is  seen  that  the  quantity  |A9(ko,  Ntb)\2  will  reach  its  maximum  growth 
value  when  f3q  =  27T (m  +  5),  where  m  is  an  integer  and  5—^0;  for  this  limit,  the  function 
sin2 (N f3q / 2) /  sin2  (f3q/ 2)  — >  N2,  i.  e.,  it  becomes  sharply  peaked  at  the  resonances  with 
increasing  N.  It  is  clear  that  this  condition  for  maximum  growth  establishes  the  selection 
rule 1'’18  for  the  photon  emission  frequency  c oq.  Indeed,  from  the  condition  f3q  =  2n m,  it 
follows  from  Eq.  (7)  in  the  radiative  long-wavelength  limit  ( qa  <C  1)  that  one  generally  has 

r  x  x  2  i1/2 


u 


q 


rnuj  b  , 


qz  qzm  —  qx 


m- 


03  B 
U3r 


(8) 


the  “Stark  ladder”  resonance  condition.  Note  that  q  is  not  collincar  with  the  direction  of 
the  applied  field  E.  then  the  integral  in  Eq.  (7)  does  not  vanish  because  of  the  periodicity  of 
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£no(k0),  and  contributes  to  the  selection  rule.  For  uq  =  ujc[1  +  (qz/Qx)2]1^2,  this  corresponds 
to  neglecting  a  small  term  ~  v±/c  as  compared  to  unity,  where  v_i_  is  the  electron  velocity 
perpendicular  to  the  y  axis.  Thus,  the  modes  that  radiate  with  the  highest  probability  cor¬ 
respond  to  the  fundamental  Bloch  frequency  and  its  harmonics.  This  quantization  condition 
is  obtained  without  requiring  any  assumptions  concerning  the  existence  of  Wannier-Stark 
energy  states. 


V.  Spontaneous  emission  probability  and  the  enhancement  factor 

The  spontaneous  emission  probability  is  evaluated  at  time  t  =  Ntb,  P*  =  Pes  (./Vtb),  by 
substituting  \Aq(k0,  Ntb)\2  from  Eq.  (6)  into  Eq.  (5).  The  sum  over  q  in  Eq.  (5)  has  been 
replaced  by  an  integral  over  q ,  taking  into  account  the  TEi0  mode  density  of  states  and 
polarization  such  that  ^  (•  •  •)  — >  (Lz/n)  f  dq( ■  ■  -)/[l  —  (qx/ q)2]1^2 ■  The  integral  can  be 
evaluated  by  using  the  property  of  the  integrand  which  contains  a  sharply  peaked,  symmetric 
function  of  q  at  q  =  qm  =  (q2  +  q2zm)1^2  [see  Eq.  (8)].  Thus,  at  every  node  defined  by  the 
resonance  conditions,  the  slowly  varying  function  of  q  in  the  integrand  can  be  replaced  by 
its  value  evaluated  at  q  =  qm,  and  then  removed  from  the  integral  over  q;  after  that,  the 
remaining  integral  can  be  evaluated  to  obtain 

q.iv^yf  K,(k0,TB)|2  (9) 

kU  [i^teAs)2]1'2 

Here  lmax  follows  from  qmax  =  lmax(wB/wc)qx,  and  determines  the  upper  limit  in  the  sum 
over  higher  Bloch  oscillation  harmonics. 

The  analysis  for  spontaneous  emission  and  radiation  characteristics  is  now  developed  by 
considering  a  SL  miniband  in  the  nearest- neighbor  tight-binding  approximation.  The  elec¬ 
tron  miniband  energy  dispersion  is  expressed  as  eno(K)  =  £no(0)  +  A  sin2 (Kya/2)  +£_l(Kj_), 
where  £„o(0)  is  the  band  edge,  A  is  the  miniband  width,  and  £_l(K^)  is  the  contribution 
from  the  perpendicular  components  of  the  miniband.  The  electron  velocity,  for  the  given 
Ky  in  the  y  direction,  is  then  given  by  vy(Ky)  =  vmaxsm(Kya),  where  vmax  =  aA/2h  is  the 
maximum  velocity  in  the  miniband. 

The  calculation  of  Pes  in  Eq.  (9)  requires  the  use  of  Hg(ko,rs)  in  Eq.  (4),  evaluated 
at  the  maximum  growth  conditions  of  Eq.  (8),  that  is  when  ftujq  =  mhus ■  In  addition, 
the  dependence  upon  q  in  Eq.  (4)  is  made  explicit  by  invoking  the  assumption  of  photon 
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long- wavelength  limit,  which  is  valid  for  all  periodic  potentials  of  interest,  even  SLs,  where 
q  <C  7T /a.  Thus,  we  find  that 

|4,1(k0,TB)|2  =  ^!fi  (10) 

Qi^b 

where  /;  =  d'OkVy('Ok)  exp(-iWk),  and  =  koya.  For  the  nearest-neighbor  tight-binding 
approximation,  one  can  find  that  //  =  —%'KVmaxd\i-  The  occurrence  of  the  Kronecker  symbol, 
<5i/ ,  allows  the  contribution  of  l  =  1  term  only,  thereby  limiting,  within  the  nearest-neighbor 
tight-binding  approximation,  the  generation  to  the  fundamental  Bloch  harmonic.  Then  it 
follows  from  Eq.  (9)  that 


PI 


2  aN 


Lx  v\ \  g1/2^c 

Ly  c 2  U)%(1  -ujI/ujI)1/2  ' 


(11) 


In  noting  that  the  SE  probability  of  Bloch  radiation  into  free  space  with  fundamental 
Bloch  frequency  is  given  by  the  expression18  Pjs  =  (27t/3 )aN  ( vm/c )2,  we  can  compare  both 
the  probabilities  analyzing  the  ratio 

,  _  p:  _ 

11  ~  pfs  klvu2b[ i  -  k2H)l1/2  ’ 

In  Fig  1,  we  show  the  enhancement  factor  r)  as  a  function  of  the  frequency  ratio  ub/xjc, 
which  determines  detuning  of  the  Bloch  frequency  with  respect  to  the  TEio  mode  cutoff 
frequency.  The  calculations  have  been  carried  out  for  a  GaAs-based  SL  (e  =  12.2)  imbedded 
into  a  rectangular  waveguide  with  horizontal  and  vertical  dimensions  Lx/ Ly  =  2.  It  is  seen 
from  the  figure  that  the  enhancement  factor  i]  increases  with  decreasing  detuning  parameter 
Ub/uJc  reaching  values  over  one  order  of  magnitude  for  the  detuning  close  to  1. 


VI.  Numerical  estimations  and  discussion 

For  numerical  estimations,  we  assume  a  GaAs-based  SL  structure  with  the  SL  lattice 
parameter  a  =  100  A,  vertical  dimension  9  /irn,  and  lateral  cross  section  18  x  1000  /mr. 
Also  we  assume  that  waveguide  is  fully  filled  with  the  semiconductor  material.  The  electron 
density  in  the  active  region  is  taken  to  be  5  x  1016  cm-3.  Taking  for  the  SL  lowest  miniband 
energy  width  A  =  20  meV,  the  maximum  group  velocity  in  the  miniband  is  estimated  as 
Vmax  =  1-6  x  107cm/s.  These  parameter  magnitudes  are  close  to  those  of  GaAs-based  SL 
structures  used  to  study  high-frequency  microwave  generation. 2,20,21  Spontaneous  emission 
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Fig.  1.  Enhancement  factor  r]  of  the  spontaneous  emission  probability 
as  a  function  of  frequency  ratio  a >b/uc  for  a  GaAs-based 
SL  structure  imbedded  into  a  rectangular  waveguide. 

of  a  photon  with  the  energy  lOmeV  corresponds  to  the  Bloch  frequency  ub  =  cub/27t  =  2.5 
THz.  The  electric  field  required  to  achieve  such  Bloch  frequency  is  E  —  Uus/ea  =10  kV/cm, 
and  results  in  the  application  of  nine  volts  across  vertical  dimension  of  the  SL  structure.23 
The  spontaneous  emission  probability  of  radiation  into  free  space  can  be  estimated  taking, 
for  example,  N  =  100  as  Pjs  =  4.3  x  10~7;  and  the  generation  energy  per  electron  hcosP/s  — 
4.3  x  10~6  meV.  Since  there  is  a  total  of  n  —  8  x  109  electrons  in  the  active  region  of  the 
SL,  the  generated  energy  achievable  is  estimated  to  be  Pfs  =  nhu>BPfs  =  34.4  eV,  which 
corresponds  to  a  power  output  generated  into  free  space  Wfs  =  ( us/N)Pfs  ~  0.14 /iW. 
Although  the  power  generated  into  free  space  is  discernibly  low  for  the  SE  of  Bloch  oscillation 
radiation,  it  follows  from  Eq.  (12)  that  SE  probabilities  and  rates  are  substantially  modified 
for  the  case  of  radiation  into  the  waveguide  mode.  We  find  from  Fig.  1  that  rj  =  20,  if 
we  take  for  the  detuning  parameter  ub/u)c  =  1-05.  Then,  using  the  obtained  value  for  the 
enhancement  factor,  we  estimate  the  power  output  generated  into  the  TEi0  waveguide  mode 
as  Wwg  ~  3  fiW. 

It  is  noted  that  a  Bloch  oscillation  SL  does  not  require  controlled  inversion  population 
between  Wannier-Stark  ladder  levels  to  get  the  desired  SE  photon  frequency;  the  desired 
frequency  is  controlled  by  the  applied  field.  Whereas  in  other  SL  light  generating  devices, 
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such  as  quantum  cascade  lasers,  a  large  inversion  population  is  required  to  provide  stimulated 
emission  with  resulting  high  threshold  current  densities  and  high  heat  dissipation.  In  this 
regard,  the  Bloch  oscillator  in  SE  offers  a  novel  option  for  operating  at  THz  frequencies, 
provided  the  power  output  can  be  enhanced  in  the  coherent  Bloch  regime.  Future  directions 
of  this  work  include  the  study  of  limiting  factors  of  electron  dephasing  due  to  scattering 
inhomogeneities  in  the  SL  to  determine  an  optimum  power  enhancement  and  the  efficiency 
of  power  extraction  from  the  cavity. 


VII.  Summary  of  the  most  important  results 

•  The  quantum  electron  dynamics  and  spontaneous  emission  of  radiation  for  a  Bloch 
electron  traversing  a  single  energy  miniband  of  a  SL  structure  in  a  constant  electric 
field  and  in  the  presence  of  a  resonant  cavity  have  been  analyzed.  The  analysis  is 
based  on  the  use  of  instantaneous  eigenstates  of  the  Bloch  Hamiltonian  good  to  all 
orders  of  the  dc  field,  and  to  first-order  perturbation  theory  in  the  quantized  radiation 
field. 

•  The  analysis  results  in  the  quantum  selection  rule  which  shows  that  the  spontaneous 
emission  into  the  waveguide  dominant  mode  is  sharply  peaked  at  frequencies  equal 
to  integral  multiples  of  the  Bloch  frequency.  This  result  has  made  use  of  no  ad  hoc 
assumptions  about  the  existence  of  Wannier-Stark  quantized  energy  levels  within  the 
band. 

•  The  spontaneous  emission  probabilities  are  substantially  enhanced  in  comparison  with 
spontaneous  emission  into  free  space  when  the  Bloch  frequency  is  tuned  by  the  field 
into  the  spectral  peak  of  the  waveguide  mode  density  of  states. 

•  The  set  of  controlling  parameters  has  been  specified  and  their  magnituds  have  been 
estimated  for  the  enhanced  spontaneous  emission  due  to  Bloch  oscillations  in  the 
terahertz  frequency  range. 

•  The  enhancement  factor  has  been  analyzed  as  a  function  of  the  Bloch  frequency  and 
the  external  electric  field.  A  theoretical  estimate  of  the  enhancement  factor  provides 
an  order  of  magnitude  enhancement  for  GaAs-based  superlattices,  and  resulted  in  a 
power  output  of  ~  3  /aW. 


9 


1  F.  Bloch,  Z.  Phys.  52,  555  (1928);  C.  Zener,  Proc.  R.  Soc.  London  A145,  523  (1934);  G.  H. 
Wannier,  Phys.  Rev.  117,  432  (1960). 

2  L.  Esaki  and  R.  Tsu,  IBM  J.  Res.  Dev.  14,  61  (1970). 

3  For  a  review  see,  for  instance,  K.  Leo,  Senricond.  Sci.  Technol.  13,  249  (1998). 

4  C.  Waschke,  H.  G.  Roskos,  R.  Schwedler,  K.  Leo,  and  H.  Kurz,  and  K.  Kohler,  Phys.  Rev.  B 
70,  3319  (1993). 

5  Y.  Shimada,  K.  Hirakawa,  M.  Odnoblioudov,  and  K.  A.  Chao,  Phys.  Rev.  Lett.  90,  046806 
(2003). 

6  K.  Unterrainer,  B.  J.  Keay,  M.  C.  Wanke,  S.  J.  Allen,  D.  Leonard,  G.  Medeiros-Ribeiro,  U. 
Bhattacharya,  and  M.  J.  W.  Rodwell,  Phys.  Rev.  Lett.  76,  2973  (1996). 

'  S.  Vinnerl,  E.  Schomburg,  S.  Brandi,  O.  Kus,  K.  F.  Renk,  M.  C.  Wanke,  S.  J.  Allen,  A.  A. 
Ignatov,  V.  Ustinov,  A.  Zhukov,  and  P.  S.  Kop’ev,  Appl.  Phys.  Lett.  77,  1259  (2000). 

8  T.  Dekorsy,  A.  Bartels,  H.  Kurz,  K.  Kohler,  R.  Hey,  and  K.  Ploog,  Phys.  Rev.  Lett.  85,  1080 

(2000). 

9  A.  W.  Ghosh,  L.  Jonsson,  and  J.  W.  Wilkins,  Phys.  Rev.  Lett.  85,  1084  (2000). 

10  N.  Marcuvitz,  Waveguide  Handbook  (Peregrinus,  London,  1993). 

11  O.  Svelto,  Pronciple  of  Lasers  (Fourth  Edition,  Plenum  Press,  New  York,  1998). 

12  L.  Mahler,  A.  Tredicucci,  R.  Kohler,  F.  Beltram,  H.  E.  Beere,  E.  H.  Linfield,  and  D.  A.  Ritchie, 
Appl.  Phys.  Lett.  87,  181101  (2005). 

13  C.  Janke,  P.  H.  Bolivar,  A.  Bartels,  H.  Kurz,  and  H.  Kiinzel,  Phys.  Rev.  B  67,  155206  (2003); 
M.  Zedler,  C.  Janke,  P.  H.  Bolivar,  H.  Kurz,  and  H.  Kiinzel,  Appl.  Phys.  Lett.  83,  4196  (2003). 

14  In  Spontaneous  Emission  and  Laser  Oscillation  in  Microcavities ,  edited  by  H.  Yokoyama  and 
K.  Ujihara  (CRC  Press,  Boca  Raton,  1995). 

15  E.  Yablonovitch,  Phys.  Rev.  Lett.  58,  2059  (1987). 

16  G.  Gallot,  S.  P.  Jamison,  R.  W.  McGowan,  and  D.  Grischkowsky,  J.  Opt.  Soc.  Am.  B  17,  851 

(2000). 

17  J.  B.  Krieger  and  G.  J.  Iafrate,  Phys.  Rev.  B  33,  5494  (1986);  G.  J.  Iafrate  and  J.  B.  Krieger, 
ibid.  40,  6144  (1989). 

18  V.  N.  Sokolov,  L.  Zhou,  G.  J.  Iafrate,  and  J.  B.  Krieger,  Phys.  Rev.  B  73,  205304  (2006). 


10 


19  J.  Singh,  Quantum  Mechanics:  Fundamentals  and  Applications  to  Technology  (A  Wiley- 
Interscience  Publication,  New  York,  1997),  p.  337. 

20  E.  Schomburg,  T.  Blomeier,  K.  Hofbeck,  J.  Grenzer,  S.  Brandi,  I.  Lingott,  A.  A.  Ignatov,  K.  F. 
Renk,  D.  G.  Pavel’ev,  Yu.  Koschurinov,  B.  Ya.  Melzer,  V.  M.  Ustinov,  S.  V.  Ivanov,  A.  Zhukov, 
and  P.  S.  Kop’ev,  Phys.  Rev.  B  58,  4035  (1998). 

21  A.  M.  Bouchard  and  M.  Luban,  Phys.  Rev.  B  47,  R6815  (1993). 

22  To  prevent  charge  domain  formation  that  may  be  possible  through  negative  differential  conduc¬ 
tance  a  superstructure  of  SLs  of  suitable  lengths  can  be  designed.  23 

23  P.  G.  Savvidis,  B.  Kolasa,  G.  Lee,  S.  J.  Allen,  Phys.  Rev.  Lett.  92,  196802  (2004). 


Accomplishments  from  this  research 

1.  V.  N.  Sokolov,  L.  Zhou,  G.  J.  Iafrate,  and  J.  B.  Krieger,  Spontaneous  emission  of  Bloch 
oscillation  radiation  from  a  single  energy  band ,  Phys.  Rev.  B  73,  205304  (2006). 

2.  V.  Sokolov,  G.  Iafrate,  and  J.  Krieger,  Spontaneous  emission  from  accelerated  Bloch 
electrons  -  Bloch  oscillation  radiation.  2006  APS  March  Meeting  (March  13-17,  2006; 
Baltimore,  MD),  Session  N  17:  Focus  Session:  Semiconductors  for  THz  and  IR  I;  Abstract 
N  17.00003. 

3.  V.  N.  Sokolov,  G.  J.  Iafrate,  and  J.  B.  Krieger,  Microcavity  enhancement  of  spontaneous 
emission  for  Bloch  oscillations ,  submitted  to  Phys.  Rev.  B  (Brief  Reports)  (2006). 


11 


Session  N17:  Focus  Session:  Semiconductors  for  THz  and  IR  I 


8:00  AM-1 0:36  AM,  Wednesday,  March  15,  2006 
Baltimore  Convention  Center  -  3 1 3 

Sponsoring  Unit:  F1AP 

Chair:  K.  K.  Choi,  Army  Research  Laboratory 

Abstract:  N1 7.00003  :  Spontaneous  emission  from  accelerated  Bloch 
electrons  --  Bloch  oscillation  radiation 

8:48  AM-9:00  AM 

Preview  Abstract 


Authors: 

Valeriy  Sokolov 
Gerald  Iafrate 

(Dept,  of  Electrical  and  Computer  Engineering,  NCSU,  Raleigh,  NC  27695-8617) 
Joseph  Krieger 

(Dept,  of  Physics,  Brooklyn  College,  CUNY,  Brooklyn,  NY  11210) 


A  theory  of  spontaneous  emission  of  radiation  for  a  Bloch  electron  traversing  a  single 
band  in  an  external  electric  field  is  presented.  The  radiation  field  is  described  by  a  free 
space  quantized  electromagnetic  field  in  the  Coulomb  gauge.  It  is  shown  that  the 
spontaneous  emission  occurs  with  frequencies  equal  to  integral  multiples  of  the  Bloch 
frequency  without  any  \textit{ad  hoc}  assumptions  concerning  the  existence  of  Wannier- 
Stark  levels.  An  explicit  expression  for  the  transition  probability  is  derived  in  first-order 
perturbation  theory  relative  to  the  radiation  field.  Although  the  output  frequency  of  the 
radiation  can  be  operationally  tuned  from  the  gigahertz  to  terahertz  spectral  range  by 
varying  the  constant  electric  field,  it  is  estimated  that  a  spontaneous  emission  power 
output  of  only  about  0.1  of  a  microwatt  is  available  using  GaAs-based  superlattices.  In 
this  regard,  it  is  noted  that  the  atomic  spontaneous  emission  probability  and  related 
transition  rates  can  be  enhanced  by  properly  tailoring  the  surrounding  electromagnetic 
environment.  Therefore,  considering  Bloch  oscillations  in  a  resonant  microcavity  to 
enhance  the  spontaneous  emission  is  a  noteworthy  alternative  for  exploring  tunable 
gigahertz  to  terahertz  radiation  sources. 


PHYSICAL  REVIEW  B  73,  205304  (2006) 


Spontaneous  emission  of  Bloch  oscillation  radiation  from  a  single  energy  band 

V.  N.  Sokolov,*  L.  Zhou,4  and  G.  J.  lafrate 

Department  of  Electrical  and  Computer  Engineering,  North  Carolina  State  University,  Raleigh,  North  Carolina  27695-8617,  USA 

J.  B.  Krieger 

Department  of  Physics,  Brooklyn  College,  CUNY,  Brooklyn,  New  York  11210,  USA 
(Received  18  December  2005;  revised  manuscript  received  8  March  2006;  published  2  May  2006) 


A  theory  for  the  spontaneous  emission  of  radiation  for  a  Bloch  electron  traversing  a  single  energy  band 
under  the  influence  of  a  constant  external  electric  field  is  presented.  The  constant  external  electric  field  is 
described  in  the  vector  potential  gauge.  The  quantum  radiation  field  is  described  by  the  free  space  quantized 
electromagnetic  field  in  the  Coulomb  gauge.  The  instantaneous  eigenstates  of  the  Bloch  Hamiltonian  are 
introduced  as  basis  states  to  analyze  the  Bloch  dynamics  to  all  orders  in  the  constant  external  electric  field.  The 
radiation  field  described  by  the  quantized  electromagnetic  field  provides  the  vacuum  fluctuations  that  are 
responsible  for  the  spontaneous  emission  of  radiation  of  the  single  electron  from  the  upper  regions  of  the 
energy  band.  It  is  shown  that  the  spontaneous  emission  occurs  with  frequencies  equal  to  integral  multiples  of 
the  Bloch  frequency  without  any  ad  hoc  assumptions  made  concerning  the  existence  of  Wannier-Stark  ladder 
levels.  An  explicit  expression  for  the  spontaneous  emission  transition  probability  is  derived  to  first  order  in  the 
quantized  radiation  field;  results  show  the  explicit  dependence  upon  the  electron  energy  band  structure,  photon 
polarization,  and  the  directionality  of  the  radiation  output.  As  an  illustration,  spontaneous  emission  probabili¬ 
ties  are  developed  and  illustrated  for  nearest-neighbor  tight-binding  and  more  realistic  superlattice  band  struc¬ 
ture  models.  For  the  GaAs-based  superlattices,  the  power  radiated  into  free  space  from  spontaneous  emission 
due  to  Bloch  oscillations  in  the  terahertz  frequency  range  is  estimated  to  be  about  one-tenth  of  a  microwatt. 

DOI:  10. 1103/PhysRevB. 73. 205304  PACS  number(s):  73.63.Hs,  72.10.Bg,  73.21.Cd,  73.50.Mx 


I.  INTRODUCTION 

Bloch  electron  dynamics  in  electric  fields  has  been  a  sub¬ 
ject  of  great  interest  dating  back  to  the  early  development  of 
solid-state  physics. 1-3  More  recently,  the  availability  of  band- 
engineered  superlattices  (SLs)  and  tailored  periodic  struc¬ 
tures  has  stimulated  further  activities4-15  in  electric  field- 
mediated  transport  and  optical  absorption  in  low-dimen¬ 
sional  SLs  and  quantum-well  (QW)  structures,  where  band- 
gaps  and  bandwidths  are  typically  several  orders  of  magni¬ 
tude  smaller  than  those  of  bulk  solids.  In  particular,  recent 
attention  has  been  focused  on  Bloch  oscillations,  the  k-space 
oscillatory  behavior  afforded  to  Bloch  electrons  when  mov¬ 
ing  in  a  periodic  energy  band  under  the  influence  of  a  con¬ 
stant  electric  field  in  the  absence  of  scattering.  For  a  biased 
semiconductor  SL,  the  Bloch  frequency1-4  for  such  an  oscil¬ 
lation  is  given  by  u>B  =  eEa/h,  where  E  is  the  applied  con¬ 
stant  field  and  a  is  the  SL  period;  it  is  evident  that  the  Bloch 
frequency  can  be  tuned  with  changing  the  external  electric 
field  so  as  to  vary  from  the  high  gigahertz  to  the  terahertz 
frequency  range,  thus  suggesting  the  applicability  of  such 
internal  oscillations  as  a  tunable  source  for  submillimeter 
wave  radiation.4 

It  is  interesting  to  note  that  transport  experiments  suggest 
only  indirect5  evidence  for  the  manifestation  of  Bloch  oscil¬ 
lations,  whereas  optical  experiments  using  semiconductor 
SLs  have  allowed  direct  observations.6,7  There  have  been 
also  numerous  discussions  relating  the  coupling  of  coherent 
Bloch  oscillations  to  other  fundamental  solid-state  excita¬ 
tions  such  as  photons,8,9  optical  phonons,10  and  plasmons.11 
In  this  paper,  in  noting  that  Bloch  oscillations  might  provide 
a  noteworthy  alternative  for  exploring  tunable  radiation 


sources  in  the  gigahertz  to  terahertz  frequency  range,  a  study 
is  undertaken  to  analyze  the  spontaneous  emission  (SE)  of 
radiation  resulting  from  a  Bloch  electron  accelerating 
through  a  single  band  under  the  influence  of  a  constant  ex¬ 
ternal  electric  field  while  subject  to  vacuum  field  fluctua¬ 
tions. 

In  the  general  framework  of  this  work,  attention  is  fo¬ 
cused  on  the  theory  for  the  SE  of  radiation  for  a  Bloch  elec¬ 
tron  traversing  a  single  energy  band  in  a  constant  external 
electric  field.  It  is  shown  that  in  a  scattering  free  environment 
and  in  fields  low  enough  to  ignore  interband  tunneling,  spon¬ 
taneous  photon  emission  occurs  as  the  Bloch  electron  inter¬ 
acts  with  the  quantum  radiation  field;  the  emission  occurs 
only  with  frequencies  equal  to  integral  multiples  of  the  Bloch 
frequency,  even  though  no  ad  hoc  assumptions  are  made 
concerning  the  existence  of  Wannier-Stark  levels.  The  tran¬ 
sition  probability  is  found  as  an  explicit  function  of  the  elec¬ 
tron  energy  band  structure,  photon  polarization,  and  direc¬ 
tionality  of  the  radiation  output. 

In  Sec.  II,  the  Hamiltonian  for  a  Bloch  electron  in  the 
quantum  electrodynamic  field  of  interest  is  developed.  The 
classical  external  electric  field  is  described  in  the  vector  po¬ 
tential  gauge,  and  the  free  space  quantized  electromagnetic 
radiation  field  is  treated  using  the  Coulomb  gauge.  In  ne¬ 
glecting  the  higher-order  quantum  field-field  interaction 
term,  it  is  shown  that  the  total  Hamiltonian  reduces  to  the 
sum  of  three  contributions,  the  Hamiltonian  for  the  Bloch 
electron  in  the  classical  external  electric  field,  the  Hamil¬ 
tonian  for  the  free  space  quantized  radiation  field,  and  the 
Hamiltonian  for  the  first-order  interaction  between  the  quan¬ 
tum  field  and  the  Bloch  electron.  In  Sec.  Ill,  the  instanta¬ 
neous  eigenstates  of  the  Bloch  Hamiltonian,  known  as  accel- 
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erated,  field-dependent  crystal  momentum  states,  are  devel¬ 
oped.15  These  instantaneous  eigenstates  are  utilized  as  basis 
states  in  describing  the  time  development,  and  in  calculating 
the  SE  transition  rates  of  the  accelerated  Bloch  electrons  un¬ 
der  the  action  of  the  perturbing  quantum  radiation  field.  In 
Sec.  IV,  the  polarization  properties  of  the  SE  resulting  from 
Sec.  Ill  are  analyzed.  Numerical  estimates  for  power  gener¬ 
ated  using  a  realistic  band  structure  are  given.  In  Sec.  V,  a 
summary  of  overall  results  is  given;  also  a  discussion  of  key 
physical  issues  relevant  to  Bloch  oscillation  dephasing  is  put 
forth,  issues  that  were  set  aside  in  the  analysis  of  this  paper, 
but  that  need  to  be  considered  in  future  work  on  this  subject. 
Also,  in  Appendix  A  we  provide  details  of  the  time-depen- 
dent  perturbation  analysis  used  to  calculate  transition  rates; 
in  Appendix  B  we  provide  the  calculation  of  matrix  elements 
utilized  in  the  perturbation  calculations  of  Appendix  A. 


II.  BLOCH  HAMILTONIAN  IN  AN  ELECTROMAGNETIC 

FIELD 


The  Hamiltonian  for  a  single  electron  in  a  periodic  crystal 
potential  subject  to  an  external  homogeneous  electric  field  of 
arbitrary  strength  in  a  time-varying  electromagnetic  field  is 

H=  t— (p  -  “A )  +  Vc(r)  +Hr.  (1) 

Here,  V,(r)  is  the  periodic  crystal  potential,  A  is  the  total 
vector  potential  consisting  of  A=AC+Ar,  and  Hr  is  the 
Hamiltonian  for  the  free  electromagnetic  field;  p  is  the  mo¬ 
mentum  operator,  r  is  the  spatial  coordinate,  m0  is  the  free 
electron  mass,  and  c  is  the  velocity  of  light  in  vacuum.  For 
the  vector  potential  terms,  Ac(?)  is  the  vector  potential  de¬ 
scribing  the  external  electric  field.  For  a  homogeneous  elec¬ 
tric  field  E(r),  turned  on  at  initial  time  t=t0,  the  vector  po¬ 
tential  is  A C{t)=—cf  \  E (t')dt' \  in  addition,  letting  F(f) 
=  <?E (f)  and  pf(f)  =  J'!(F(f')c/f',  it  then  follows  that  pc(f) 
=  -(e/c)Af(r).  The  vector  potential  term,  A,,  corresponding 
to  the  free  quantized  radiation  field,  is  given  as 


A,.= 


■(flq je“ 


qr 


+  alf 


/q*r 


), 


(2) 


where  eqj  is  a  unit  polarization  vector  for  the  radiation  mode 
with  wave  vector  q  and  polarization  j,  and  eq  /- eq  /(  =  <5;// 
with  j  J'  =  1,2.  aq  .  and  «q  /  are  the  creation  and  annihilation 
boson  operators  of  the  quantum  radiation  field,  V  is  the  vol¬ 
ume  of  the  system,  and  </  =  |q|,  the  magnitude  of  wave  vector 
q.  A,,  satisfies  the  Coulomb  gauge,  that  is,  V-A,.=0,  which 
means,  from  Eq.  (2),  that 


Sq,;  •  q  =  0,  (3) 

for  each  polarization  direction.  It  also  follows  from  the  Cou¬ 
lomb  gauge  that  [p,A,.]  =  0.  The  Hamiltonian  for  the  free 
quantized  radiation  field  is 

Hr=^fi(oqa]qjaqj,  (4) 

q  J 

where  u>q=cq,  the  free  space  photon  dispersion. 


In  substituting  the  vector  potential  A=At,+A).  into  Eq. 
(1),  and  regarding  the  term  containing  A,,  as  a  perturbation,16 
we  can  reduce  the  total  Hamiltonian  of  Eq.  (1)  to 

H  =  H0  +  Hr  +  Hh  (5) 

where 

H0(t)  =  “  [p  +  pe(f)]2  +  Cc(r) ,  (6) 

2/71q 

and 


H,(t)  =  -  —  A,.  •  [p  +  p(,(r)] ,  (7) 

m0c 

where  Hr  is  given  by  Eq.  (4).  Equation  (5)  follows  from  Eq. 
(1)  in  the  limit  of  a  low  intensity  radiation  field,  where  we 
have  dropped  the  term  of  order  A2.17  In  Eq.  (5),  the  first  two 
terms  represent  the  Hamiltonian,  //0(f),  for  a  single  electron 
in  a  periodic  crystal  potential  interacting  with  a  homoge¬ 
neous  electric  field,  and  the  Hamiltonian,  Hn  for  the  free 
electromagnetic  field;  the  term  Hj{t)  in  Eq.  (5),  and  specified 
in  Eq.  (7),  couples  both  subsystems  H0  and  Hn  and  causes 
transitions  between  the  accelerated  Bloch  electron  states 
through  photon  absorption  and  emission.  In  the  next  section, 
starting  with  the  total  Hamiltonian  of  Eq.  (5),  use  is  made  of 
first-order  time-dependent  perturbation  theory  to  calculate 
SE  transition  probabilities  between  states  of  H0+Hr  while 
regarding  Hj{t)  as  a  perturbation. 

III.  QUANTUM  DYNAMICS  BASED  ON 
INSTANTANEOUS  EIGENSTATES 

A.  Crystal  Hamiltonian  and  instantaneous  eigenstates 

In  this  analysis,  the  dynamical  properties  are  considered 
for  the  situation  in  which  the  electron  is  confined  to  a  single 
band  “n0”  of  a  periodic  crystal  with  energy  e„o(K);  the  ef¬ 
fects  of  interband  coupling15  and  electron  intraband  scatter¬ 
ing  are  ignored.  Therefore,  the  quantum  dynamics  is  de¬ 
scribed  by  the  time-dependent  Schrodinger  equation, 

ihj[\%,o(t))=H\vno(t)y,  (8) 

here,  H  is  given  by  Eq.  (5),  and  |\P,„  (f))  is  sought  in  terms  of 
the  complete  set  of  eigenstates  based  on  the  instantaneous 
eigenstates  of  H0  and  the  eigenstates  of  Hn  subject  to  an 
initial  Bloch  momentum  state  and  initial  photon  field  occu¬ 
pancy.  The  instantaneous  eigenstates  of  H0  satisfy  the 
equation15 

(  ^“[P  +  PfM]2  +  ye(r)  j  <A,ok(f)  =  e„o[k(r)](/f„ok(f),  (9) 

where 

eiKr 

<An0k(r)(r,f)  =  —^unoHt){r,t)\  (10) 

here,  un  k(()(r,f)  is  the  periodic  part  of  the  Bloch  function, 
k(t)  =  K  +  (l/^,)/JoF(r')^t'=K  +  p(..(?)/^,  where  the  values  of 
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K  are  determined  by  the  periodic  boundary  conditions  of  the 
periodic  crystal  of  volume  O.  The  basis  states  for  Hr  are 
given  by  the  well-known  free  photon  field  equation, 

Hr\{nqJ})  =  2  h(oqnqJ\{nqj}) ,  (11) 

q  J 

where  |{nq  /})  is  a  simple  product  of  all  possible  combina¬ 
tions  of  photon  number  states,  nqp  with  a  given  wave  vector 
q  and  polarization  sq  /. 

Thus,  the  solution  to  |'vI,„(j(f))  of  Eq.  (8)  can  be  represented 
in  terms  of  eigenstates  of  basis  states  | t/r„  k((),{nqj}) 
-|i/,H0k(/))|{,7qj})  °f  the  unperturbed  Hamiltonian  H0+Hr  as 


It  is  clear  from  above  that  |A^  *[k0(r)  ,t]\2  is  the  probability 
for  the  emission  of  a  photon  with  wave  vector  q  and  polar¬ 
ization  eq j  at  time  t.  Yet,  the  total  emission  probability  as 
calculated  in  first-order  perturbation  theory  results  in  a 
weighted  sum  with  respect  to  (nq  -+1)  when  summed  over 
all  photon  wave  number  q  and  polarization,  as  noted  in  Eq. 
(13).  Therefore,  it  follows  that  if  no  photons  are  initially 
present  in  the  quantum  radiation  field,  that  is  nq  ;  =  0  for  all 
modes,  then  the  emission  process  results  in  the  well-known 
SE  probability, 

^W  =  2|<](k0,f)|2.  (15) 

q  J 


l^„0(?))  =  E  E  A{nqj}(k,?)l^ok»,{"qj}) 

k  { nqj> 

Xexp(_  J '  [en0[k(r')]  +  2  h(oqnqJ\dt' 

(12) 

where  the  summation  over  k  is  carried  out  over  the  entire 
Brillouin  zone  (BZ),  and  {nqj}  is  specified  over  all  possible 
combinations  of  photon  occupation  number  with  correspond¬ 
ing  photon  wave  vectors  q  and  polarization  eq  ;. 

The  appropriate  time-dependent  equations  of  motion  for 
the  }(k,f)  coefficients  expressed  in  Eq.  (12)  are  given  in 
Appendix  A;  these  equations  relate  the  time  dependence  of 
Airl  i(k ,  t)  to  the  basis-dependent  matrix  elements  of  H, 
through  a  self-consistent  set  of  equations.  In  applying  the 
usual  first-order  time-dependent  perturbation  theory  method¬ 
ology  to  this  set  of  equations,  consistent  with  the  initial  con¬ 
ditions  for  the  k  state  of  the  Bloch  electron  in  the  energy 
band  “n0”  and  the  initial  state  of  radiation  field,  a  first-order 

perturbation  theory  result,  namely,  a!"  ,(k,r)=Alfl)  ,(k,f) 

,  ,  '"qjt  t"q j> 

+A,  ,(k,t),  is  obtained,  where  the  superscripts  “a”  and  “e” 

\nq,js 

stand  for  photon  absorption  and  emission,  respectively.  The 
explicit  expressions  for  both  Aj''1  j(k,r)  and  A^’1  Tk,f)  are 
given  in  Appendix  A,  in  Eqs.  (A6)  and  (A7). 

In  the  case  of  photon  emission,  the  total  emission  prob¬ 
ability  can  be  written  as 

^w=EKJ+i)K:»(k0,r)i2,  (i3) 

q  J 

where 

D  C1 

Aq)(k0,0  =  -7=  dt'  v[k0(f')  -  q]  •  eq  j 
Vg-Lo 

-yj  K[Mt)] 

J  f0 

-e„0[k0(fi)  -qj-fcwjdfj,  (14) 

with  k0(f)=K0+pc(t)/^.  Here,  K0  and  nqj  are  the  initial  val¬ 
ues  of  the  wave  vector  K  and  photon  number  nq  j,  and  D 
=—i'\2'nalV,  where  a=e2/(hc)  is  the  fine  structure  constant. 


The  additional  term,  proportional  to  nqj,  in  Eq.  (13)  is  also 
well  known  to  correspond  to  induced  emission  probability. 

The  terms  associated  with  photon  absorption  in  Eq.  (A6) 
can  be  considered  in  a  similar  way;  here  we  find,  for  the 
absorption  probability,  that 

^W  =  2«qJKj(k od)\\  (16) 

q  J 

where 
<j(  k„,r) 

D  f 

=  —  dt’  v[k0(r')  +  q]  •  eqJ 

v qJ>0 

Xexp(-  j  J '  [e„0[k0(fi)]  -  e^tk^fj)  +  q]  +  ha)q\dt 

(17) 

The  quantity  Pa(t)  in  Eq.  (16)  is  proportional  to  nqj,  and 
therefore  vanishes  if  the  radiation  field  is  initially  in  the 
vacuum  state. 


B.  Properties  of  transition  probabilities  and  selection  rules  for 
spontaneous  emission  of  a  single  photon 

The  specific  analysis  now  focuses  on  one-photon  sponta¬ 
neous  emission  from  a  Bloch  band  “n0”  coupled  to  a  radia¬ 
tion  field  that  initially  has  no  photons  present.  This  means 
that  the  probability  amplitude  in  the  wave  function  of  Eq. 
(12)  satisfies  the  initial  condition  A{„  j(k,t0)={<5„  ,o}^k,k0 
at  time  t=t0  when  the  electric  held  is  turned  on.  With  this 
assumed  initial  condition,  the  probability  amplitude  for  emis¬ 
sion,  Aj](k0,f),  is  now  evaluated  from  Eq.  (14). 

In  evaluating  A^'.(k0,f),  it  is  assumed  that  the  external  dc 
held,  E,  is  along  the  z  axis  (also  the  growth  direction  of  a 
SL);  then,  it  follows  that  k0z(t)  =  K0z+F(t-t0)/h,  and  k0±(f) 
=  K0±  =  const,  where  k0±(f)  is  the  component  of  wave  vector 
k0(f)  perpendicular  to  the  z  axis.  In  taking  advantage  of  the 
periodic  properties  of  the  terms  in  Eq.  (14),  A^.(k0,f)  is 
evaluated  in  clocked  integral  multiples  of  the  Bloch  period, 
so  that  t=NrB,  where  tb  =  2it/ a>B,  the  time  to  traverse  one 
period  of  the  BZ.  Then,  the  integral  in  Eq.  (14)  over  time  can 
be  replaced  by  an  integral  over  k0z  through  the  substitution 
dr=(h/ F)dk0z.  It  then  follows  from  the  periodic  property  of 
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the  energy  band  in  the  BZ,  and  therefore  the  electron  veloc¬ 
ity  as  well,  that  the  probability  amplitude,  at  integral  mul¬ 
tiples  of  the  Bloch  period,  can  be  expressed  as15 


A«(k0,MrB)  = 


1  ~  exp(-  iN/3q) 
1  -  exp(-  ipq) 


(18) 


Here,  the  parameter  fiq  is  given  by 

[  rK oz+Gz 

Pq  =  -  I  [s„0(k0)  -  e„o(k0  -  q)  -  h(oq]dk0z,  (19) 

J  K0z 

and  A^.(k0,  tb),  the  probability  amplitude  over  a  single 
Bloch  period,  is 


hD  fKo?+G? 

=  ~r  dk0z  v(k0  -  q)  •  Sq j 

F<qJKQz 

Xexp(  -  J  J  [fi/i0(ko)  -  e«0(ko  -  q)  -  h(»q]dk'0z  j , 

(20) 


where  G,=2tt/ci,  the  z  component  of  the  SL  reciprocal- 
lattice  vector.  It  then  follows  from  Eq.  (18)  that 


|AW(k0,lVr,)|2  =  ^flA-Cko,.)!2.  (2D 


From  Eq.  (21),  it  is  seen  that  |A^j(k0, A^T/j))2  will  reach  its 
maximum  growth  value  when  (3q=2Tr(m  +  S),  where  in  is 
an  integer  and  5— >0;  for  this  limit,  the  function 
sin2(N(3q/2)/sm2(/3q/2)^>N2.  It  is  clear  that  this  condition 
for  maximum  growth  establishes  the  selection  rule 15  for  the 
photon  emission  frequency.  Indeed,  from  the  condition  j3q 
=  2  mu,  it  follows  from  Eq.  (19)  that 


j  rKoz+Gz 

(oq  =  ma)B+—  [e„(k0) 

n°zJ  K0z 


en0(k0  -  q)]^0z- 


(22) 


V1 


COn 


FIG.  1.  Dependence  of  the  relative  probability  spectral  density 
y{coq)  of  SE  on  the  normalized  photon  frequency  coql coB. 


r  Koz+cz 

v±(K±)  =  (l/Gj  vjk0)dk0z,  (24) 

Jk0z 

with  v±(k0)  =  (l/^)VK^e„o(K)|ko.  Specihcally  for  the  tight- 
binding  energy  band  approximation,18  one  finds  that 
Vj_(K_l)  =  vj_(Kj_).  However,  the  second  term  on  the  right- 
hand  side  of  Eq.  (23)  is  negligibly  small  for  nonrelativistic 
Bloch  electron  velocities;19  thus  one  generally  has 


u>„  —  me oB,  q  —  qm  =  m — ,  (25) 

c 

the  “Stark  ladder ”  resonance  condition.15 

Figure  1  illustrates  the  resonance  frequency  behavior  of 
the  relative  probability  spectral  density  r]((i)q)  =  |A|'|(k0, 
Ntb)\2/\A^  j(k0,  tb)\2  for  N=4.  With  increasing  N,  the  r/icOy) 
becomes  sharply  peaked  at  the  resonances  because  the  func¬ 
tion  v(jr)  =  sin2(AiTrjr)/[A?  sin2 (tot)]  behaves  like  a  delta  func¬ 
tion  at  each  N  as  0.  Thus,  the  modes  that  radiate  with  the 
highest  probability  correspond  to  the  fundamental  Bloch  fre¬ 
quency  and  its  harmonics;  this  quantization  condition  is  ob¬ 
tained  without  requiring  any  assumptions  concerning  the  ex¬ 
istence  of  Wannier-Stark  energy  states. 


In  noting  the  contribution  from  the  integral  in  Eq.  (22),  note 
that  if  q  is  collinear  with  the  direction  of  the  applied  field, 
i.e.,  in  the  z  direction,  then  the  integral  in  Eq.  (22)  vanishes 
because  of  the  periodicity  of  s„o(k0).  If,  however,  q  is  not 
collinear  with  the  direction  of  the  applied  field,  then  the  in¬ 
tegral  in  Eq.  (22)  does  not  vanish,  and  contributes  to  the 
selection  rule.  In  particular,  in  considering  the  radiative  long- 
wavelength  limit,  where  sn  (k0-q)  can  be  expanded  in  a 
Taylor  series  about  small  q,  the  expression  for  u>q  in  Eq.  (22) 
results  in 


coq  =  mcoB  +  q±  •  v±(K±) .  (23) 


C.  Total  spontaneous  emission  probability 

The  total  SE  probability  is  evaluated  from  Eq.  (15).  The 
SE  probability  is  evaluated  at  time  t=NrB  by  substituting 
|Aq)(k0 ,^tb)|2,  already  obtained  in  Eq.  (21),  into  Eq.  (15)  to 
obtain 


;  P].(Nth) 


V 


(2tt) 


Cl  max 

;  dqq ■ 

Jo 


.sin2(A/y2) 
sin2(/3  /2)  J0 


r47r 

dn%  l< 

Jo  i 


(26) 


Here,  q  is  the  perpendicular  component  of  q  relative  to  the 
applied  electric  field  direction  and  the  average  transverse  ve¬ 
locity,  v±,  is  given  by 


Here,  the  sum  over  q  in  Eq.  (15)  has  been  replaced  by  an 
integral  over  q  for  a  single  polarization,  such  that  Sq(-  •  •) 
^[V/(2rr)3]/ dq  q2 f  dfl(-  ■■),  where  <ffl  =  sin  Q dO dip  is  the 


205304-4 


SPONTANEOUS  EMISSION  OF  BLOCH  OSCILLATION... 


element  of  a  solid  angle  subtended  by  q,  and  6 ,  <p  are  the 
polar  angles.  In  treating  the  integral  of  Eq.  (26),  it  is  ob¬ 
served  in  the  integral  that  tb) |2  is  a  slowly  vary¬ 

ing  function  of  q\  whereas,  it  is  noted  from  Fig.  1,  and  pre¬ 
vious  discussions  therein,  that  the  term  sin2(A'/i(//2)/ 
sin2(/3f//2)  is  a  sharply  peaked,  symmetric  function  of  q  at  q 
values  of  qm=mcoB/ c,  where  m  is  an  integer,  as  noted  in  Eq. 
(25).  Thus,  at  every  node  defined  by  the  resonance  condi¬ 
tions  /3q  =  2mn,  with  u>ql u>B=m  and  qm=rna>B/ c,  the  slowly 
varying  function  of  q  in  the  integrand  can  be  replaced  by  its 
value  evaluated  at  q  =  qm,  and  then  removed  from  the  integral 
over  q:  the  remaining  term  in  the  integrand  can  be  evaluated 
by  letting  q  =  (a>B/2irc)f3q  at  each  node  so  that 

f  sin2(M/y2)  _  r  r/£,sin2(jV/y2)  a>g 

J  q  sin2(/39/2)  7 TcJ_n  2  sin2(/?9/2)  c 


Thus,  Eq.  (26)  becomes 


Pi 


<»b  V 
c  (27 r)3 


dn2\A^0,rB)\\  (27) 

j 


where  Zmax  follows  from  <7max=ZmaxcoB/c,  and  determines  the 
upper  limit  in  the  sum  over  higher  Bloch  oscillation  harmon¬ 
ics. 

The  calculation  of  Pse  in  Eq.  (27)  now  requires  the  use  of 
A^.(k0,TB)  in  Eq.  (20),  evaluated  at  the  maximum  growth 
conditions  of  Eq.  (23),  that  is  when  hcoq=mh(oB 
+  /iq  i  •  v  |  (K  | ).  In  addition,  the  dependence  upon  q  in  Eq. 
(20)  is  made  explicit  by  invoking  the  assumption  of  a  photon 
long-wavelength  limit;  this  assumption  is  valid  for  all  peri¬ 
odic  potentials  of  interest,  even  SLs,  where  q<Tr/ a.  Thus,  in 
letting  e„0(k0-q)==e  (k0)-q-Vkoe„o(ko)  and  v(k0-q) 
—  v(k0)-q-  Vk(v(k0),  A^.(k0,  tb)  in  Eq.  (20)  becomes,  in  the 
long-wavelength  limit, 

Aq  j(k0,  tb ) 

fin  rKoz+Gz 

=  —e‘^K0z>  dk0 ,  v(k„)  •  eq  , 

Fiq  J  k0z 

(28) 


with  the  phase  factors  in  the  exponent  defined  by 


i p(x)  =  a 


(mx-qz 


e„0W  \ 

huB  / 


(29a) 


x(koz,  K(h) 


a 

u>B 


(k0z-  K0z)q± 


(29b) 


The  results  of  Eqs.  (28)  and  (29)  are  valid  for  an  arbitrary 
band  structure.  But  for  the  tight-binding  band  structure 
model  cases  now  considered,  especially  since  they  are  appro¬ 
priate  for  SLs,  it  is  noted  that18  Vj_=v_l  in  Eqs.  (29b),  in 
which  case  x{kl):,K{)?)  =  {).  Therefore,  for  the  tight-binding 
model,  Eq.  (28)  can  be  simplified  by  introducing  dimension¬ 
less  variables  and  parameters. 
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enn(^0;) 

& k  =  k0za ,  =  qza,  eB(&k)  =  — - .  (30) 

na>B 

Then  the  integral  in  Eq.  (28)  becomes 

<j( k0,rB)  =  -  ^e'W<Dm.(#9),  (31) 

uBiq 

where 

‘M#,)  =  ^f  dVkTJVk)e-‘m\  (32) 

the  mth  Fourier  component  of  the  integrand  rq  ;(0t,) 
=V8qj<"V»*',*(  Thus,  it  follows  from  Eq.  (31)  that 

|Aq,j(k0,  rB)\2  =  (33) 

wBq 

where  D  is  defined  after  Eq.  (14).  Substituting  Eq.  (33)  into 
Eq.  (27),  one  obtains 

tlj  ^max  r Att 

K  =  —^qi\  dtiZ  |$//#9|)|2.  (34) 

Mb0  1=1  J  o  j 

In  analyzing  the  polarization  properties  of  the  SE,  it  is 
convenient  to  separate  the  contribution  to  with  po¬ 

larization  parallel  and  perpendicular  to  the  applied  electric 
field,  i.e.,  relative  to  the  z  axis,  so  that  <bm/(  iTf/)  =  <\>mj7(  i)q) 
+  where 

=  (e, j)z^~  [  d$k  vz(\)e^W-^ 

(35) 

and 

^mjJV  =  '  2tj-J  d®k v±Weil*‘>sd*khm*ki- 

(36) 

Here  the  symbol  J_  stands  for  the  vector  component  perpen¬ 
dicular  to  the  applied  electric  field.  These  terms  determine 
angular  dependences  of  the  intensity  of  SE.  The  next  step  is 
to  evaluate  the  integrals  in  Eqs.  (35)  and  (36),  which  requires 
a  specific  model  for  the  energy  band  structure.  This  will  be 
addressed  in  the  next  section. 

IV.  ANALYSIS  OF  SPONTANEOUS  EMISSION  FOR 
SUPERLATTICE  AND  TIGHT-BINDING  MODELS 

The  analysis  for  spontaneous  emission  and  radiation  char¬ 
acteristics  is  now  developed  by  utilizing  the  special  case  of  a 
SL  miniband  with  a  growth  direction  defined  along  the  z 
axis.  The  electron  energy  band  dispersion  relation  is  ex¬ 
pressed  as 

00 

e  (K)  =  e  (0)  +  2  A,  sin2  ^  +  Ae(K±) ,  (37) 

/=i  1 

where  e„o(0)  is  the  band  edge.  A/  is  the  width  of  the  Zth 
miniband  harmonic  of  the  SL,  and  Ae(K±)  is  the  contribu- 
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tion  from  the  perpendicular  components  of  the  band.  The 
corresponding  velocity,  for  the  given  Kz  in  the  z  direction,  is 
then  given  by  v  Z{KZ)  =  (l  /  h)\_deno(Kz)  /  dKz]=^=iv  /  sin (Kzla), 
where  U/=A//a/2^,  the  maximum  velocity  associated  with 
the  /th  miniband  of  bandwidth,  A/.  The  energy  band  disper¬ 
sion  of  Eq.  (37)  in  the  SL  direction  generally  includes  long 
range  coupling  over  the  neighboring  QWs  with  a  relative 
strength  measured  by  the  specific  value  of  the  ratio  A/+1/A, 
<  1 ,  which  is  dependent  upon  the  extent  of  wave  function 
overlap.  For  the  well-known  case  of  nearest-neighbor  tight- 
binding  (NNTB)  energy  dispersion,  only  A!  is  considered 
significant,  so  that  next  nearest  neighbor  and  longer  range 
QW  wave  function  overlaps  are  assumed  to  be  negligibly 
small.  Also  for  the  energy  band  model  under  consideration  in 
Eq.  (37),  it  is  noted  that  v±(k0)  is  independent  of  k(]z,  so  that 
v±  =  v±  (Ref.  18)  from  Eq.  (24),  and  for  subsequent  use  in 
Eq.  (36). 

A.  Direction  of  spontaneous  emission  radiation  output 
relative  to  the  electric  field 


FIG.  3.  Geometry  of  radiation  for  longitudinal  SE  (d=0). 


=  (y  •  SqjMm.O-  (38) 

Thus,  for  L  polarization,  only  the  l=m  term  survives,  which 
is  in  resonance  with  the  mth  Bloch  harmonic  oscillation.  In 
contrast,  the  model  with  NNTB,  which  is  frequently  used  to 
approximate  electron  transport  in  narrow  energy  bands, 
would  allow  the  generation  of  /  =  1  only,  thereby  limiting  the 
propagation  to  the  fundamental  Bloch  harmonic. 

Substituting  <t,„y,,(0)  from  Eq.  (38)  into  Eq.  (33),  the  L 
polarization  probability  amplitude  becomes 

T  2 

|A^i(k0,rB)|2  =  (39) 

V  qu~B 

where  RT(qk)  is  introduced  as  an  angular  form  factor.  The 
emission  probability  density  is  independent  of  the  angle  i p, 
the  angle  the  q  makes  with  the  x  axis  in  the  xy  plane  (see 
Fig.  2).  Thus,  the  form  factor  Rjiip),  denoted  by  R-r( & 
=  tt/2  ,  cp)  is  unity  on  the  right-hand  side  of  Eq.  (39).  Thus, 
the  SE  probability  is 


uN  ^  , 

- - 2  qpl  Sfl 

/=1 


max 

Y3  V I 

1=1  c 


(40) 


First,  it  is  useful  to  analyze  the  particular  cases  of  SE 
probability  for  the  q  direction  of  propagation,  where  q  is 
perpendicular  to  E,  the  transverse  radiation  direction,  and 
then  where  q  is  parallel  to  E,  the  longitudinal  radiation  di¬ 
rection.  Finally,  this  will  be  followed  by  the  case  where  q  is 
at  an  arbitrary  direction  with  respect  to  the  electric  field. 

1.  Radiation  in  transverse  direction 

For  the  case  of  the  photon  wave  vector,  q,  perpendicular 
to  the  applied  field,  E,  then  as  noted  from  Fig.  2,  qz= 0  and 
9=|q±|;  <4so>  the  two  independent  polarization  directions, 
eq ;,  can  be  chosen  parallel  to  the  field  direction,  with  longi¬ 
tudinal  polarization  [referred  to  as  L  polarization,  with  j=  1 
and  (eq  l).=  l,  (eq  l)±  =  0],  and  perpendicular  to  the  field  di¬ 
rection,  with  transverse  polarization  [referred  to  as  T  polar¬ 
ization,  with  7=2  and  (eq2):=0,  |(eq 2)xl  =  1]-  Since  q.= 0,  it 
follows  from  Eq.  (30)  that  =  0  so  that  Eqs.  (35)  and  (36) 
become  analytically  manageable  for  evaluation,  and  one  ob¬ 
tains 

oo 

^my,z(0)  —  —  ~  (^q V lfy,m  ~  ~  ~  (^qj/z^mj 
1  1=1  L 


where  account  has  been  taken  for  photon  wave  vector  ip 
=  la>B/c  and  the  small  solid  angle  SD,  =  2tt  SO  for  0e[(7t 
-S0)/2,(tt+S0)/2\,  S0<Z  1. 

The  emission  of  photons  with  T  polarization  gives  a  zero 
contribution  into  Pse,  since  the  condition  (m=l=0,j=2)  in 
Eq.  (38)  results  in  zero  wave  vector  qt= 0=0. 

2.  Radiation  in  longitudinal  direction 

For  the  case  of  the  photon  wave  vector,  q,  parallel  to  the 
applied  field,  E,  while  the  polarization  eqj  is  perpendicular 
to  the  direction  of  the  field,  then  consider  q  =  \qz\,  q  =  0  with 
(f5q/).  =  ()  and  |(eq/)±|=  1 .  The  two  independent  polariza¬ 
tions,  eq  /,  are  chosen,  as  noted  in  Fig.  3.  In  this  case,  since 
#qA0,  the  term  ,  (iff  in  Eq.  (36)  can  be  calculated 
analytically  if  we  omit  the  higher  harmonic  contribution  (/ 
>  1)  in  the  energy  band  dispersion  of  Eq.  (37).  Then  Eq.  (36) 
becomes 

=  (v  •  eqj)  j  exp(—  iirm/2) ,  (41) 

where  Jm(x)  is  the  Bessel  function  of  the  first  kind. 
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Thus,  from  Eq.  (41)  the  longitudinal  component  of  SE 
will  be  nonvanishing  for  finite  values  of  Vj_;  the  electron  just 
moves  along  the  in-plane  direction  at  a  constant  Vj_,  which  is 
fixed  by  the  initial  quasimomentum.  Substituting  Eq.  (41) 
into  Eq.  (33),  the  probability  amplitude  becomes 


Kj(ko,TB)|2  =  « 


(2tt)3  1  (v  •  eqJ)i  r2  l  #„A, 


V  q  col  \2hb 


Rl(<p)- 


(42) 


Then,  from  Eq.  (34),  the  SE  probability  becomes 


ocN 


Lily  j  VT  9 

K= — 

wbc  ;=i 


2  (  1l 

2hc 


m. 


(43) 


where  <5I7  =  27r[l -cos(<56?)]~  tt(SO)2,  the  small  solid  angle 
about  the  z  axis,  where  S0<  1.  The  photon  polarization  per¬ 
pendicular  to  the  velocity  v ,  does  not  contribute  to  the  SE. 
RL(cp),  the  emission  form  factor  in  Eq.  (42),  is  isotropic  in 
the  xy  plane  of  Fig.  3;  therefore  Rl(ip)  =  Rl(0=(),<p)=  I .  By 
utilizing  the  fact  that  ql=la>B/c  in  the  argument  of  the  Bessel 
function  in  Eq.  (43),  it  follows  that  ftqAl/(2fuoB)  =  lvl/c 
<S  1.  Thus,  the  asymptotic  form  for  Jt(x)  —  (x/2)1  /\'(l  +  I ).  for 
1  is  assumed,  where  T (x)  is  the  gamma  function,  and 
v 2xxa!2h\  therefore,  Pse  in  Eq.  (43)  reduces  to 

^max  ;2/+l  /  \  2/ 

m  <44) 


In  particular,  for  the  fundamental  Bloch  harmonic  (/=  1),  Eq. 
(44)  becomes 


K 


2 

(v±/c)2  SCI. 


(45) 


3.  Radiation  in  arbitrary  direction 

For  the  case  of  photon  wave  vector,  q,  at  an  arbitrary 
angle  with  respect  to  the  electric  field,  as  noted  in  Fig.  4(a), 
the  angular  dependence  of  the  SE  probability  can  be  eluci¬ 
dated  from  Eqs.  (35)  and  (36).  Since  the  integrals  of  Eqs. 
(35)  and  (36)  can  only  be  evaluated  numerically  when  the 
full  energy  dispersion  of  Eq.  (37)  is  used,  use  is  made  of  the 
NNTB  approximation  in  Eq.  (37)  by  keeping  the  term  with 
/=  1  only  in  the  dispersion  relation.  In  this  case,  the  longitu¬ 
dinal  component,  $my-z(#?),  in  Eq.  (35),  is  found  to  be,  after 
integration 

m  (Ofj 

=  (eqjjj - Jm(v  ^-/wB)exp(-  imnl 2)  .  (46) 

clz 

In  noting  again,  as  in  subsection  IV  A  2,  that  the  argument  of 
the  Bessel  function  is  small  so  that  an  asymptotic  limit  may 
be  imposed,  it  follows  that  Eq.  (46)  reduces  to  <I,,„;  ,(i^/) 
=  (eq  j)zv1m(vlqJ  exp(-i7rm/2)/(2"1w!). 

The  transverse  component,  |  ( dq) ,  in  Eq.  (36)  can  be 
calculated  in  a  similar  fashion.  It  is  found  that 
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FIG.  4.  (a)  Geometry  of  radiation  for  an  arbitrary  wave  vector  q. 
(b)-(c)  Form  factor  Rj(6,(p)  of  the  SE  probability;  (b)  fundamental 
Bloch  frequency  (/=  1);  (c)  second  Bloch  harmonic  (1=2). 

\®mj(®q)\2  =  \(Eqj)zm—  +  (v-EqjM  J2m{v xqJ wB) . 

(47) 

For  wave  vectors  q  with  a  z  component  projected  over  qz 
=qZm  =  (m(oB  cos  6)1  c,  the  second  term  in  the  square  bracket 
is  small  compared  to  the  first  term  with  respect  to  the  param¬ 
eter  y=|vjJ/c,  except  in  a  narrow  angular  range  defined  by 
S9—  y.  This  correction  can  be  taken  into  account  through  the 
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function  gmj(0,  <p)  =  1  +gz(v’eqj)x/[wwB(eqj)J,  which  is 
equal  to  1  if  the  latter  term  is  ignored.  Then  we  obtain  from 
Eqs.  (34)  and  (47), 


K 


dil 
cos2  0 


,21  Vl 
J )  ]  / —  COS 


0  2  (eq j)2zgfj(e,<p) . 

'  7=1 


(48) 


For  iq/  sin  01  c  1 ,  this  equation  takes  the  form 

?max  ,2/+l  /  \  2/  r4l7 

''■'“"Siwlc)  J„  <49) 


where  the  form  factor 


fl/(0-<p)  =  cos2(/  *’(0)2  (®qj)^L(0>^)  (50) 


2=1 


determines  the  angular  dependence  of  the  SE  probability  for 
the  /th  Bloch  harmonic.  In  particular,  at  zero  transverse 
velocity  [g/j(0, <p)  =  l],  we  obtain  the  form  factor,  taking 
into  account  the  fundamental  (/=1)  and  the  second  (1=2) 
Bloch  frequency  harmonics  /?i(0)  =  sin2(0)  and  R2(  0) 
=  (1 /4)sin2(2$),  respectively.  For  nonzero  transverse  veloc¬ 
ity,  the  form  factor  depends  on  both  polar  angles  (0,(p),  with 
the  dependence  on  <p  being  weak  according  to  the  function 
gij(0,ip).  Figures  4(b)  and  4(c)  illustrate  the  anisotropy  of 
the  form  factor  corresponding  to  the  fundamental  (/  =  1)  and 
the  second  (1  =  2)  Bloch  harmonic.  Finally,  integrating  over 
all  polar  angles,  we  find 


(51) 


B.  Spontaneous  emission  discussion;  Power  estimate 

The  spontaneous  emission  probability  has  been  obtained 
and  analyzed  for  a  variety  of  directional  outputs.  In  general, 
it  is  clear  from  Eq.  (34)  and  (47)  that  the  average  probability 
per  unit  time,  resulting  from  first-order  time-dependent  per¬ 
turbation  theory,  is  Pse(t=NTB)/(NTB)~  u>B;  in  contrast,  in  a 
one-electron  atom,  with  a  two-level  model,  the  resulting  fre¬ 
quency  dependence12  for  the  rate  of  SE  ~<u3g,  where  coeg  is 
the  characteristic  frequency  corresponding  to  the  electron 
transition  between  the  excited  state  and  the  ground  state. 
Thus,  the  SE  probability  of  Bloch  radiation  has  a  compara¬ 
tively  weaker  frequency  dependence  than  the  corresponding 
two-level  atomic  SE  probability. 

Finally,  a  numerical  estimate  is  performed  for  purposes  of 
evaluating  the  power  generated  during  the  spontaneous  emis¬ 
sion  process.  In  this  regard,  a  semiconductor  SL  of  one  mi¬ 
cron  length  is  considered,  consisting  of  100  periods  with  SL 
lattice  parameter  a  =  100  A,  and  with  lateral  cross  section  S 
=  1 00  X  100  /i m2.  The  SL  is  assumed  to  have  an  electron 
density  of  5  X  1017  cm-3,  which  corresponds  to  a  total  num¬ 
ber  of  electrons  in  the  SL  active  region  of  72  =  5  X  109.  For  the 
SL  miniband  structure,  in  the  NNTB  approximation,  it  is 
assumed  that  the  lowest  miniband  energy  width,  Ah  is 


20  meV,  so  that  the  maximum  group  velocity  in  the  mini¬ 
band  [/=  1;  Eq.  (37)]  is  V\=aA1/(2h)=1.6X  107  cm/s;  these 
parameter  magnitudes  resemble  those  of  GaAs-based  SLs 
used  to  study  high-frequency  microwave  generation.4,20'21 

Spontaneous  emission  is  considered  for  a  photon  energy 
given  by  hu>c=fiu>B=\Q  meV,  which  corresponds  to  the  fun¬ 
damental  Bloch  frequency  vb=o>b/(2tt)=2.5  THz.  The  cor¬ 
responding  electric  field  required  to  achieve  this  Bloch  fre¬ 
quency  is  E=fiu>B/(ea)=  10  kV/cm,  and  results  in  the 
application  of  one  volt  across  the  SL.  It  then  follows  from 
Eq.  (51)  that  the  SE  probability  for  N=  100  is  estimated22  to 
be  Pse(N=  100)  =4.3  X  10-7.  From  this  SE  probability,  an  es¬ 
timate  of  the  generation  energy  per  electron  becomes 
ha>BPse= 4.3X  1 0  6  meV,  and  since  there  are  a  total  of  n  =  5 
X  109  electrons  in  the  active  region  of  the  SL,  the  generated 
energy  achievable  is  estimated  to  be  P =nhu>BP[  =  2\ .7  eV, 
which  corresponds  to  an  approximate  power  generation  of 
W=(vb/N)P  —  0A  /2.W.  In  this  power  estimate,  it  is  noted 
that  the  power  generated  in  W  is  proportional  to  coB.  and  thus 
E2,  as  expected  from  power  considerations.  Although  the 
power  generated  is  discernibly  low  for  SE  of  Bloch  oscilla¬ 
tion  radiation  into  free  space,  it  is  noted  that  SE  probabilities 
and  rates  can  be  modified  by  tailoring  the  surrounding  reso¬ 
nant  electromagnetic  environment. 12-14  It  would  therefore  be 
interesting  to  consider  the  prospect  of  enhancing  superlattice 
spontaneous  emission  of  Bloch  oscillation  radiation  through 
superlattice-resonant  cavity  interaction  with  the  cavity  tuned 
to  the  Bloch  frequency,  while,  at  the  same  time,  admitting 
scattering  processes  as  an  offsetting  suppressant;  this  is  the 
subject  of  a  future  investigation. 

V.  SUMMARY  AND  DISCUSSION 

A  theory  for  the  spontaneous  emission  of  radiation  for  a 
Bloch  electron  traversing  a  single  band  in  a  uniform  electric 
field  and  in  a  scattering-free  environment  has  been  devel¬ 
oped.  The  analysis  results  in  a  selection  rule  that  shows  that 
spontaneous  emission  is  sharply  peaked  at  frequencies  equal 
to  integral  multiples  of  the  Bloch  frequency,  a  result  that 
made  no  ad  hoc  assumptions  concerning  the  existence  of 
Wannier-Stark  quantized  energy  levels  within  the  band,  but 
that  comes  directly  from  the  use  of  electric  field-dependent 
instantaneous  eigenstates  of  the  Bloch  Hamiltonian.  Sponta¬ 
neous  emission  probabilities  were  derived  and  analyzed  in 
terms  of  band  structure  and  polarization  radiation  properties. 
A  theoretical  estimate  for  the  intensity  of  the  spontaneous 
emission  for  GaAs-based  superlattices  resulted  in  a  power 
output  into  free  space  of  one-tenth  of  a  microwatt. 

The  theoretical  analysis  has  implicitly  assumed  a  single 
electron  model  for  a  Bloch  electron  in  a  uniform  electric 
field  of  a  perfect  superlattice  crystal  ignoring  scattering  ef¬ 
fects.  Such  assumptions  provide  the  optimal  conditions  for 
coherent  Bloch  oscillations  and,  therefore,  reflect  the  condi¬ 
tions  for  the  maximum  achievable  power  output,  albeit,  one- 
tenth  of  a  microwatt  into  free  space,  as  calculated  in  Sec.  IV. 
But  in  recent  years,  many  studies  on  coherent  Bloch  oscilla¬ 
tions  have  pointed  to  the  importance  of  additional  scattering 
effects  such  as  carrier-carrier  scattering,25  LO-phonon  scat¬ 
tering,26  and  alloy  disorder  (interface  roughness)  scattering27 
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due  to  superlattice  compositional  doping,  all  effects  that 
strongly  influence  the  dephasing  of  coherent  Bloch  oscilla¬ 
tions,  and  therefore  have  a  significant  influence  on  the  mag¬ 
nitude  of  the  spontaneous  emission  output.  The  dephasing 
effects  result  in  a  broadening  of  the  peaks  of  the  spontaneous 
emission  probability  function  that  determines  the  selection 
rule  for  transitions,  in  our  case  cocl=mcoB,  and  also  dampen 
the  output  THz  radiation  amplitude.  Specifically,  interface 
roughness  scattering27  has  been  identified  as  a  dominant 
mechanism  for  dephasing  in  GaAs/AlGaAs  superlattices 
over  a  wide  temperature  range;  competitive  dephasing  mech¬ 
anisms  have  been  reported  due  to  LO-phonon  emission  in  a 
miniband  wider  than  the  LO-phonon  energy  (36  meV  for 
GaAs-based  superlattices),  and  also  from  carrier-carrier  scat¬ 
tering  effects.25,26  In  all  cases,  the  estimated  linewidth  for 
broadening  from  such  effects  ranges  as  —(1.25-2.5)  meV, 
which  corresponds  to  dephasing  times  in  a  range  of 
—  (1.0-0. 5)  ps;25  27  this  broadening  in  energy  is  nontrivial  as 
it  corresponds  to  about  20%  of  the  Bloch  frequency. 

It  is  also  noted  that  the  external  electric  field  in  our  analy¬ 
sis  is  assumed  to  be  homogeneous  across  the  active  region  of 
the  superlattice.  However,  in  a  realistic  device  analysis,  as 
the  electron  density  is  increased  to  achieve  maximum  effi¬ 
ciency  for  spontaneous  emission,  the  Bloch  oscillating  super¬ 
lattice  can  develop  inhomogeneous  fields.  In  this  case,  the 
developed  theory  can  be  extended  by  treating  the  Bloch  elec¬ 
tron  dynamics  in  inhomogeneous  electric  fields  of  arbitrary 
strength  and  time  dependence.15,28  In  other  situations,  do¬ 
main  formation  may  be  possible  through  negative  differential 
conductance.  In  this  case,  by  superlattice  design,  a  super¬ 
structure  can  be  tailored  as  a  stack  of  superlattices  of  appro¬ 
priate  lengths  to  prevent  domain  formation.14 

Last,  it  is  noted  that  a  Bloch  oscillation  superlattice  does 
not  require  controlled  inversion  population  between  Wan- 
nier-Stark  ladder  levels  to  get  the  desired  spontaneous  emis¬ 
sion  photon  frequency;  the  desired  frequency  is  controlled  by 
the  applied  field.  Whereas  in  other  superlattice  light  generat¬ 
ing  devices,  such  as  quantum  cascade  lasers,  a  large  inver¬ 
sion  population  is  required  to  provide  stimulated  emission 
with  resulting  high  threshold  current  densities  and  high  heat 
dissipation.  In  this  regard,  the  Bloch  oscillator  in  spontane¬ 
ous  emission  offers  a  novel  option  for  operating  at  THz  fre¬ 
quencies,  provided  the  power  output  can  be  enhanced  in  the 
coherent  Bloch  regime. 

In  closing,  the  power  emitted  into  free  space  from  the 
spontaneous  emission  of  Bloch  oscillation  radiation  is  dis- 
cernibly  small.  It  is  our  intention  in  a  future  effort  to  consider 
enhancing  superlattice  spontaneous  emission  of  Bloch  oscil¬ 
lation  radiation  through  superlattice-cavity  inetaraction  with 
the  cavity  tuned  to  the  Bloch  frequency,  while,  at  the  same 
time,  including  the  offsetting  effects  of  dephasing  inhomoge¬ 
neities  mentioned  earlier,  so  as  to  more  realistically  evaluate 
the  optimal  magnitudes  of  power  output  from  the  spontane¬ 
ous  emission  of  Bloch  oscillation  radiation. 
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APPENDIX  A:  RESULTS  OF  PERTURBATION  THEORY 
APPROACH 

The  equation  for  the  probability  amplitudes  Ain  }(k ,  t)  in 
Eq.  (12)  is  given  by 


dt  in  , ,  r  ,  , 

k 

X */'„0k(/)|tf/l  40k'(r),{«qj}> 

Xexp(-i[,0. 


e„0[k'(?')]-e)!o[k(r')] 


+  2  {n'qJ  ~  nqJ)h(oq 
q  J 


(Al) 


We  assume  that  at  initial  time,  t0,  the  system  is  in  one  of  the 
eigenstates  of  Hamiltonian  H0+Hr  with  wave  function 
liA»0K0’{)1q  .}),  corresponding  to  the  Bloch  electron  in  a  single 
band  “n0”  with  the  wave  vector  K0,  i.e.,  ij/„  k0 

=  ( 1  / nl/2)e'K»''u„nKn.  and  with  the  initial  distribution  of  pho¬ 
ton  numbers  in  the  radiation  field  |{/;6  ■}).  Substituting 

A{nll/k ,  ,}(k ’ ;)  ,}(k , f)  +  -  •  int0  eT  (A1)’  and 

taking  into  account  the  initial  condition  A{„  i(k,f0) 
=  { Sn  „o  }4(f),k0(r)>  one  obtains  to  the  zeroth  and  first  order 
in  Hb  defined  in  Eq.  (7),  for  a!0)  ,(k,f)  and  Al1*  ,(k,r),  re- 

\wqjj  \nq  ,j! 

spectively, 

A^.}(k,t)  =  {4q.,„Oj}4(f),ko(r),  (A2) 


where  k0(f)  =  K0+pf(f)/ft  and  {<5.  ,,o  }  =  EL  ,<5,,  ,.o  .  Matrix 

u'  q j  q J'  q j 

elements  for  perturbation  operator  Hj  (Appendix  B)  are 
given  by 


({Wq,/K  iA«0k(/)l^/l  *A;iQk0(f)>{wq.y}) 


27ra 

V 


2  ^fVksJk(t)] 


q  j 


X  W  ;{4qT"l})40-k-q' 


+  VV/'  +  1(4,,.„  n°,  .,+U  {4  ,«°  ))4„,k+q'), 

'  11  w  q  J  q  j  q -J  q ,/  0  H 


(A4) 


where  a=e2l(hc)  is  the  fine  structure  constant.  Thus,  we 
obtain 
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<A5> 

where  the  first  term  on  the  right-hand  side  of  Eq.  (A5)  is  due 
to  absorption  of  a  photon  from  the  radiation  field  by  the 
Bloch  electron 


q  j 


x4 

Kl 


v£  r  d, 
>-k-  v  L 


v[k(f')] 


^  J‘o 

Xexp Mk°(0] 


■  e«0[k(f")]  +  ho)q,]df  ) , 


(A6) 


and  the  second  term  represents  emission  of  a  photon  by  the 
Bloch  electron  to  the  radiation  field 


<^.<>=02  v-w  + 1(4,. 


q  j 


x4 


f  dt’ 


v[k(r')] 


v<?  -o 
Xexp^-  ;  |  {e„o[k0(r")] 

-  e„Jk(r")]  -  hw  ,}df  ) . 


(A7) 


Here  v[k(r)]  =  (l//f)VKe„o(K)|k(f)  is  the  instantaneous 
velocity  of  the  Bloch  electron  confined  to  the  band 
“«o”  D=-i^2Tra/V,  and  (Sn  ,,  „o  ±1;{<5„  ,„o  }) 

q  J  ’  q'j  qj’  q j 

=  4,  .,,n°,  ,±inqj^q»j/5n  -0  . 

q  J  q'j  H  J  ^  J  <U  q j 

Using  these  results  and  Eq.  (A5),  we  can  easily  calcu¬ 
late  |a|*  j(k,f)|2.  When  calculating  the  double  sum 
2q  j  2q  j-  (...),  we  make  use  of  the  relationships  resulting 
from  properties  of  the  Kronecker  delta  symbol.  Then 
the  crossing  terms  vanish 


A<fl)*  (k,f)A|e)  ,(k,f) 

\nqjr  '  iwq,/T  ’  7 


=A)  '  i(k,f)A; ;  Tk,f)=0.  Thus,  we  obtain 

iwqjj 

|Ag/M|2  =  l<.}(M|2  +  \A$Jk,t)\2,  (A8) 


where 


q'j' 

•  P  dr'  v[k(t')] 
V<7  -Ln 


Xexp 


-  7  f  {s„0[k0(f")]  -  e„o[k(f")] 

nj>o 


and 


+  fuOq’jdt" 


2tt 


(A9) 


I4v}(k’f)l2= 2  (v,t  + 


X  4„.k+q' 


vv  K 


v[k(r')] 


Xexp 


"  *1 K 


[k0(r")] 


■  E„0[k(r")]  -  h<aqi]df 


(A  10) 


APPENDIX  B:  MATRIX  ELEMENTS  OF  PERTURBATION 
HAMILTONIAN 

Matrix  elements  of  the  perturbation  Hamiltonian  are  cal¬ 
culated  by  substituting  Hj  from  Eq.  (7)  and  ^llQk(t)  °f  Eq. 
(10),  while  using  the  properties  of  the  photon  annihilation 
and  creation  operators,  («qj-l|aqj|nq;)  =  (nqj)1/2  and  (nqj 

3qjlnq-^“^,q J" 


+  l|a'  :\n„  j)=(naj+ 1)1/2.  The  result  is 


(lVq./} 5  'A„(jkU,|W/l  *An,jk ' f fj ? { ^ q j'} ) 


e  2irfic  eq'J' 


moc  V  E  V?' 


x  V'V/'(4v,v,  ,-i;{4  v )) 

\  q  q  v  q',/  qj  qj 

x  J  ^«0k(P  +  Pc)e,q'rfA„0k'  dr 


V^q'  jr  !(4nt  ,+l’i4  j,nr  .}) 

q  J  q  j  q'j'  qj 


/ 


x  k(p  +  Pc)e  'q  r'A„(1k'  dr 


(Bl) 


In  using  iA,,0k(p  of  Eq.  (10),  the  first  integral  is  equal  to  zero, 
except  when  K'-K  +  q  =  G,  where  G  is  the  vector  of  the 
reciprocal  lattice  (Here  we  omit,  for  brevity,  the  prime  on  the 
q  vector).  For  values  of  K  and  q  in  the  first  BZ,  we  can  take 
G  =  0.  Then,  it  follows  that  K'=K-q,  and  since  k(f)=K 
+  p c{t)lh,  then  we  get  k'  =k-q.  In  this  case,  we  have,  for  the 
first  integral  in  (Bl), 


££'SJ"  ^0k(P  +  P< 

J 


Vqr<A„nk'  dr 


=  4',k-qJ  <y  °W„0 k)  (p  +  pt)«  °%k-q*0, 

(B2) 

where  the  integration  over  dr()  is  carried  out  over  the  primi- 
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tive  cell  volume.  In  the  long-wavelength  limit  (qa<6  1),  /j^, 
can  be  approximated  as 

4tk'  =  4',k-q  j  (e/Kr°M„0k)*(P  +  Pc)e‘Kr°un0k  *0 
=  4',k-q  J  l/f«0k(P  +  Pc) l40k  dr .  (B3) 

In  noting  that  the  well-known  momentum24  expectation 
value  for  Bloch  states  is 

J  (e'k(')rM„ok)*p(e'k(,)ru„ok)r/r  =  m0v[k(f)],  (B4) 

where  v(k)  =  (l/ft)Vk£„o(k),  it  then  follows,  for  k(r)  =  K 
+  p C(t)/H  in  Eq.  (B4),  that 


'/fI0k(P  +  Pc)  <A>0k  dr  =  w0v(k) ,  (B5) 

where  </fJI()k  are  the  instantaneous  eigenstates  of  Eq.  (10). 
Thus,  Eqs.  (B3)-(B5)  become 

7ktk'  =  4\k-q'«ov(k) .  (B6) 

Similarly,  for  the  second  integral  in  Eq.  (Bl),  one  finds  that 

7k~k'  =  J  <A«0k(P  +  Pr)e~‘<ir'/\k’  dr  =  4',k+q™ov(k) . 

(B7) 

Then,  using  (Bl),  (B6),  and  (B7),  the  matrix  elements  for  the 
perturbing  Hamiltonian,  Hh  are  established  for  use  in  Eq. 
(A4). 
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